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We analyse the quantum behaviour of the “Little Sibling” of the Big Rip singularity (LSBR) 
[I[. The quantisation is carried within the geometrodynamical approach given by the Wheeler- 
DeWitt (WDW) equation. The classical model is based on a Friedmann-Lemaitre-Robertson- 
Walker Universe filled by a perfect fluid that can be mapped to a scalar held with phantom character. 
We analyse the WDW equation in two setups. In the first step, we consider the scale factor as the 
single degree of freedom, which from a classical perspective parametrises both the geometry and 
the matter content given by the perfect fluid. We then solve the WDW equation within a WKB 
approximation, for two factor ordering choices. On the second approach, we consider the WDW 
equation with two degrees of freedom: the scale factor and a scalar field. We solve the WDW 
equation, with the Laplace-Beltrami factor-ordering, using a Born-Oppenheimer approximation. 
In both approaches, we impose the DeWitt (DW) condition as a potential criterion for singularity 
avoidance. We conclude that in all the cases analysed the DW condition can be verified, which 
might be an indication that the LSBR can be avoided or smoothed in the quantum approach. 
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I. INTRODUCTION 

From the early thinkers that devoted their reasoning to the nature of motion and the related conceptions of space 
and time, to modern theories of gravitation, geometry has always been an essential aspect of the epistemological 
bridge between observations and theory. General relativity (GR) opened the door for a dynamical role of space-time 
geometry in physics. 

GR was the main theoretical foundation for the development of relativistic astrophysics and modern cosmology. 
Accordingly, the evolution of black hole research inspired the formulation of the singularity theorems of Hawking and 
Penrose [H (together with the cosmic censorship conjecture, excluding the possibility of naked singularities i). These 
efforts contemplated past and future singularities. The first, associated with the Big Bang, were supposed to obey the 
Weyl curvature hypothesis, following thermodynamical motivations Q, whereas future singularities in that time were 
of two types: astrophysical (gravitational collapse of a star into a black hole) and cosmological (big crunch scenarios). 
On the other hand, the remarkable evolution of research in cosmology, both theoretically and observationally, led 
to the most accepted paradigm of (AGDM) cosmology (with inflation), which still contains the unresolved puzzles 
of dark matter (DM) and dark energy (DE) M- The empirical data supportiM a present accelerated expansion 
(see for example 0), motivated the re-appearance of the cosmological constant [j-d,[Tlj. Nevertheless, the absence 
for a corresponding physical theory (or an inconsistency with the quantum field predictions for the vacuum energy) 
remains. Moreover, a constant dark energy density gives the so called coincidence problem, i.e. the same order of 
magnitude at present time for the DM and DE densities. These issues inspired alternative descriptions. Dynamical 
dark energy models based on scalar fields were widely investigated and the research is still active. 

From the point of view of physical ontology, the ideal study would follow from a fundamental theory describing 
the form of the equation of state for a DE fluid, based on the understanding of its origin and constituents. In fact, 
different cosmological evolutions can be explored depending on the nature of the energy-momentum tensor and on 
the equation of state (satisfying or transgressing some energy conditions). Until the present date, the nature of dark 
energy remains unclear, therefore cosmologists have been somewhat limited to testing different equations of state (and 
different potentials for the hypothetical scalar fields), and compute the resulting cosmic dynamics. In this context, 
dynamical DE models based on scalar helds have many applications. Some models try to unify DM and DE as in 
the case of the generalised Ghaplygin gas [l^ - [l^ and also early-time and late-time acceleration, see for example [l^ . 
Other models assume an interaction between the DM and DE components (see [l^ and 0 for recent observational 
constraints). But most importantly, was the discovery that dynamical dark energy models known as phantom energy 
(which might have a negative kinetic term) could give rise to a new type of singularity, the Big Rip |18l - l26l | , in which 
the scale factor, the Hubble parameter and its cosmic time derivative diverge in a finite cosmic time. It was the 
beginning of the study of late-time or future singularities, or more generally, of cosmic singularities related to dark 
energy Other types of singularities were discovered such as. Sudden singularity [T7I-0. Big Freeze [sol - l^ . type IV 
[30l - l^ . Little Rip [381 - 0 and the recently investigated Little Sibling of the Big Rip (LSBR) [I|. 

Future cosmic singularities or abrupt events are fascinating new areas of cosmology. Gan these events be absent 
when the classical theory is quantised? Some works have ap plied the methods of quantum cosmology to the classical 
models leading to these DE related singularities 0. 143 - 0 . In fact, dynamical dark energy models based on scalar 
fields (quintessence or phantom) provide a very interesting and suitable scenario to explore the quantum gravity 
challenge. Accordingly, in quantum cosmology the restriction to an isotropic and homogeneous universe, simplihes 
substantially the general theory. The approach is done in the line of quantum geometrodynamics where a canonical 
quantisation of gravity (metric functions and conjugate momenta) is performed |54l. l55|. The main purpose is to solve 
the Wheeler-DeWitt (WDW) equation and apply appropriate boundary conditions, in order to get the wave function 
of the Universe [54 - 10 . In it was shown that after solving the WDW equation with phantom dark energy with 
a corresponding exponential potential, the quantum effects dominate the region of the classical Big Rip singularity, 
where classically the scale factor, the Hubble rate and it cosmic time derivative blow up at a finite future cosmic time. 
These authors found wave packets solutions that follow the classical trajectory and disperse in the genuinely quantum 
region. These quantum effects occur at large scales and since the solutions are regular the Big Rip singularity is 
considered to be effectively avoided in the quantum analysis. The quantum cosmology of the classical Big Brake 
singularity was analysed in [45j |. This is a type of Sudden future singularity where the Hubble rate reaches zero and 
the deceleration approaches infinity leading to an abrupt brake of the expansioiQ. It was found that under reasonable 
assumptions, the DeWitt (DW) criterium is satisfied, i.e. the wave function vanishes in the region of the classical 
singularity. The Big Demurrage, another kind of Sudden singularity, and the Big Freeze, the later can be seen as a Big 


^ In fact, a Sudden singularity occurs at a finite scale factor where the Hubble rate is finite but its cosmic time derivative diverges. 
Depending on the parameters of the equation of state, this singularity can happen in the future as a Big Brake l45l or in the past. 
This second case, it what we named the Big Demarrage, it results from a phantom fluid leading to a cosmic expansion that starts with 
infinite acceleration [^ . 
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Rip happening at a finite scale factor, were also investigated in the quantum approach based on the WDW equation 
1^ . These singularities can result from a dark energy fluid with an appropriate generalised Chaplygin gas equation 
of state M- In the last mentioned cases (including type IV singularity), the DW criterium for singularity avoidance 
is satisfied, pointing to a possible avoidance of the singularities An essential result from all these works is 

that singularity avoidance in quantum cosmology necessarily predict quantum effects at scales much larger than the 
Planck length. 

The heart of quantum cosmology is to apply the quantum theory to the Universe as a whole. This contrasts 
with the approach in which there is a classical theory with quantum effects (corrections) in certain phases of the 
evolution. Therefore, quantum cosmology must be based on a quantum theory of gravity. Three major promising 
approaches are String theory. Loop Quantum Gravity and Quantum Geometrodynamics, where the last two are 
different developments of Canonical Quantum Gravity [^ . Independently of the correct theory of quantum 
gravity, the quantum cosmology based on the WDW equation should remain valid at least on energies below the 
Planck scale (if not on all scales) [s^. The two main challenges of quantum cosmology consist in i) finding a 
description of the dynamical evolution and ii) in finding the quantum state of the universe with the appropriate initial 
or boundary conditions. These issues are interrelated although we do not have a robust theory for the boundary 
conditions or the initial quantum state Nevertheless, the WDW equation is crucial for the topic of boundary 

conditions in quantum cosmology. A relevant feature of the WDW equation is the fact that, for usual scalar fields, it 
is locally hyperbolic (taking the form of a wave equation) such that it has a well posed initial value problem. This 
fact is related to the indefinite sign of the kinetic term, which is directly linked to the attractive nature of gravity 
for usual matter [s^. The presence of a phantom field changes the structure of this equation, for example, if there 
is “phantom dominance”, the WDW equation becomes elliptic (or parabolic), which changes the imposition of the 
boundary conditions. This is of extreme relevance since while the solutions of hyperbolic equations are “wave-like”, a 
perturbation of the initial (or boundary) data of an elliptic (or parabolic) equation is felt “at once” by essentially all 
points in the domain. Another important aspect is related to the problem of time: the WDW equation is independent 
of an external time parameter. On the other hand, the origin of the arrow of time can in principle be related to the 
structure of this equation [H, . 

In summary, the quantum cosmology of dark energy models leading to classical singularities, is showing very clearly 
that quantum effects are predicted at scales much larger than the Planck length and that it is possible to find solutions 
to the WDW equation that could effectively avoid the classical singularity. Many open questions remain, such as the 
appropriate boundary conditions, the classical-quantum correspondence, the problem of time and the interpretation 
of the wave function. 

In this work, we apply the methods of quantum cosmology to the LSBR [l[ . The quantisation of the classical model 
presented in [l| will be carried within the geometrodynamical approach given by WDW equation. We consider the 
case where the LSBR is induced by a perfect fluid, in which the scale factor is the single degree of freedom, and 
also the case where the LSBR is induced by a scalar field with phantom character; therefore in a cosmological model 
with two degrees of freedom. In the canonical quantisation we have [a, iTa] ^ 0, where tTo is the canonically conjugate 
momentum of the scale factor a, therefore in principle one can choose different “factor-ordering” in the derivation 
of the WDW equation [sj. Accordingly, for the case with a perfect fluid we will solve the WDW equation within a 
WKB approximation, for two different factor ordering choices. For the phantom scalar field case, we solve the WDW 
equation within the Laplace-Beltrami factor-ordering [sj, with two degrees of freedom: the scale factor and a scalar 
field, which classically describe the geometry and the matter content respectively. The solutions are obtained using a 
Born-Oppenheimer (BO) approximation. In all the approaches, we impose the DW condition and find that it can be 
verified which might be an indication that the classical LSBR event can be avoided in the quantum description. As 
will be pointed in the conclusions, the DW condition on the wave function might not be sufficient to guarantee that 
all the relevant quantities which diverge in the classical model become finite in the quantum description. One should 
compute the corresponding expectation values and/or probability amplitudes of the physically relevant quantities to 
have a complete analysis. In this work, we simply check if the DW condition can be verified. The calculation of the 
expectation values and the probability amplitudes requires further investigation since it is intimately linked to various 
open questions in quantum cosmology regarding the Hilbert space structure of the solutions and the classical-quantum 
correspondence. 

The structure of the article is as follows: in section II, we review the basic aspects of the classical model analysed 
in [I| and introduce a phantom scalar field which could induce a classical LSBR. Section III includes the quantum 
analysis for the case of a perfect fluid where the scale factor is the only degree of freedom; i.e. the LSBR is induced 
by a perfect fluid which can be fully determined through the scale factor which equally defines the geometry. We 
review the WDW equation, find its solutions using the WKB approximation and verify the DW condition using two 
different factor ordering. In section IV we solve the WDW equation for the phantom scalar field case, using a BO 
approximation and evaluate the asymptotic behaviour of the corresponding solutions to see if the DW condition can be 
satisfied. Finally, in section V we present the conclusions. Three appendixes are also included: Appendix A contains 
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some mathematical details of the WKB approximation used in section III and also the validity of this method. Some 
detailed calculations related to the Parabolic Cylinder functions appearing in section IV are relegated to appendix B. 
Finally, in appendix C we present the conditions for the validity of the BO approximation used in section IV. 


II. THE LITTLE SIBLING OF THE BIG RIP EVENT: REVIEW 


The classical model leading to the future event we are interested in was exposed in Ref. [I|, nevertheless, we present 
here the basic features of the classical description relevant for the present work. We consider the Einstein equations 
compatible with the cosmological principle, therefore the Universe is assumed to be homogeneous and isotropic, i.e. 
it is well described by a Friedmann-Lemaitre-Roberston-Walker (FLRW) space-time metric and the matter content 
is represented by a perfect fluid. The Friedmann equation governing the evolution of the scale factor is 


SttG k 

— Ptot - 


( 2 . 1 ) 


where k = —1, 0,1 for a hyperbolic, flat or spherical spatial geometry respectively, H is the Hubble parameter and ptot 
includes different forms of matter-energy, such as matter (baryonic-l-dark-matter), radiation and some form of unknown 
DE. We disregard the interaction between these fluids, therefore all obey the conservation equation V = 0, where 
is the energy-momentum tensor for each component and stands for the covariant derivative. Since this implies 
that for pressureless matter Pmat ~ while for radiation we have Pi-ad then the asymptotic behaviour of 

the scale factor will be governed by the dark energy component pde- Note that for the late-time evolution, we can 
neglect the curvature term. In fact, for a dominating form of dark energy which by definition is transgressing the 
strong energy condition (SEC: p -I- 3p > 0), the DE density, pde will always dominate over the curvature term. The 
cosmological constant transgresses the SEC condition, but there are many other possibilities coming from a whole 
range of models. For these type of fluids, the Universe faces an accelerated expansion, according to the Raychaudhuri 
equation 


a 

a 


47rG 


(p + 3p). 


( 2 . 2 ) 


where p is the pressure and the dot stands for derivative with respect to the cosmic time. As in Ref. [ij we assume 
an equation of state that slightly deviates from that of a cosmological constant (in the following sections we drop the 
“de” label) 


P = -P-^- (2-3) 

We will consider A to be positive and very small such that the model can mimic closely the ACDM in the past and 
during the present [l|. For this choice, the conservation equation for p(a) gives a logarithmic dependence with the 
scale factor [I| 


p{a) = A + Aln(j-y (2.4) 

where ag is an integration constant. In [I|, ag was regarded as the present value of the scale factor such that p(ao) = A 
mimics a cosmological constant at present. In this work, although we are using the same notation as in [l|, we are 
interested in the late-time cosmological evolution, therefore we will consider the solution of the Friedmann equation 
for the late-time expansion during a DE dominated epoch. For the fluid obeying the equation of state (ESI), the scale 
factor diverges in the inhnite future as well as the Hubble rate and the four-dimensional Ricci curvature, but the 
derivative of the Hubble rate does not. This future event is what we refer as “Little Sibling” of the Big Rip singularity 

Q- 

The dynamics of the previously described perfect fluid can be represented by a scalar held with a phantom character; 

i.e. 




= -y 





where p^ = p and = p. Using the equation of state (12.31) we get 


d(j) ^ rK 1 

da V 3 aH (a) ’ 


V{a)=A+^ 


l-b61n 



(2.5) 


( 2 . 6 ) 
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FIG. 1. Plot of the scalar field versus the logarithmic scale factor a = ln(a/ao) as given by the expression (12.811 . The little 
sibling of the big rip is located at large values of a. 


and considering the (asymptotic) Friedmann equation, 


= 


SttG 


-P, 


(2.7) 


with p defined in Eq. (1^ . we get two sets of general solutions for the classical trajectory in configuration space: 

1 


(j){a) = ± 


IA 




In 


where (/>i is an integration constant fixed as 


(j)i = (j){ao) =F 


IA 


+ 01 ) 


A' 


( 2 . 8 ) 


(2.9) 


The evolution of the scalar field in terms of a = In is depicted in Fig. [T] Finally, equation (12.81) can be inverted 
and through Eq. (EH) we get a uni-parametric family of quadratic potentials V (0) 


^(0) — "T + 27rAG (0 — 0i)^ , 
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( 2 . 10 ) 


which is shown in Eig. [51 


III. QUANTUM STUDY WITH A PERFECT FLUID 
A. Wheeler—DeWitt equation - review 

The WDW equation can be introduced in the canonical quantisation of gravity in the so called quantum ge- 
ometrodynamics approach, we next summarise how it can be obtained under the assumptions of homogeneity and 
isotropy. From the Einstein-Hilbert gravitational action functional, the Einstein (Euler-Lagrange) equations can be 
obtained and used to deduce the Hamiltonian equations describing the dynamical evolution. It turns out that this is 
a constrained dynamical system, in which the constraints correspond to the invariance of the model with respect to 
space-time diffeoinorphisms. There are two types of constraints: the momenta constraints (space diffeomorphisms) 
and the Hamiltonian constraint (time diffeomorphisms) 1^. 

For cosmological applications it is usual to consider the principle of homogeneity and isotropy; i.e, FLRW space-time 
metric. The Hamiltonian is derived from the Lagrangian, from which the dynamical equations can be obtained. In 
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FIG. 2. Plot of the potential (12.1011 . The small vertical line close to the origin corresponds to the initial value </>!. The abrupt 
little sibling of the big rip happens at large values as shown by the black arrows. The circles indicated schematically the 
occurrence of the little sibling of the big rip. 


fact, the classical Hamiltonian constraint 'H = 0 corresponds to the Friedmann equation. The passage to the quantum 
description is done by promoting all the degrees of freedom and canonically conjugate momenta to operators acting 
on some Hilbert space. The usual approach is to consider the so called minisuperspace, the space of all 3-dimensional 
totally symmetric (FLRW) metrics and all matter-energy configurations. It is analogous to the configuration space 
in the classical analysis. The wave function of the universe (in the Schrodinger representation) will live on this space 
and each point on it represents a FLRW universe with a certain value of the scale factor and certain well defined 
values for all the matter-energy degrees of freedom. The WDW equation governs precisely the distribution of the wave 
function in this (quantum) configuration space. It corresponds to the Hamiltonian constraint (time diffeomorphism 
invariance) previously mentioned, where the Hamiltonian operator acting on the wave function gives zero. Therefore, 
time is absent in this quantum description and the wave function in the Born interpretation is expected to provide 
a stationary probability (amplitude) distribution in minisuperspace [s^. As mentioned in the introduction, in the 
canonical quantisation procedure we have to take into account that there are different “factor-orderings” in obtaining 
the WDW equation. 

The action for gravity with a FLRW symmetry can always be written in the following way [H, [s^, 

Stt f 

S = Sg + Sm, Sg = — I [-a‘^a +ka] dt, (3.1) 

where Sm corresponds to the matter part which can be represented by a perfect fluid dependent on the scale factor (one 
single degree of freedom) or include intrinsic matter-energy degrees of freedom like a scalar field, (j). The corresponding 
WDW equation can be expressed as 

nip = 0, (3.2) 

where T-L is the quantum Hamiltonian computed from the above action. 

A reasonable model is one in which the decoherence of the general superposition of quantum states is absent in 
those regions of minisuperspace corresponding to a classically behaving universe . This requirement is analogous to 
the correspondence principle in usual quantum mechanics where the classical regime is re-obtained for large quantum 
numbers. Here, coherence presupposes the possibility for confined wave packet solutions, effectively reproducing the 
classical trajectory in configuration space. However, in the vicinity of classical singularities, where the quantum 
cosmology can have important non-classical effects, decoherence might be present and the quantum results can differ 
from the classical description [i^. Therefore, for any FLRW model with singularities (or abrupt events such as the 
LSBR), one should check if the distribution of the wave function in minisuperspace is such that these (classical) cosmic 
events can be effectively avoided with the quantum approach. 

In this section, we study the quantum behaviour of the LSBR induced classically by the perfect fluid with equation 
of state ()2.3I) . We investigate the solutions of the WDW equation and apply the DW criterium for effective singularity 
avoidance, namely, the condition that the wave function vanishes in the region corresponding to the classical LSBR. 
We will consider also two different factor ordering choices. 
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B. Wheeler—DeWitt equation with a perfect fluid - WKB approximation 


We start by considering the simplest case within the WDW equation where the matter content is given by a perfect 
fluid with a given equation of state, therefore whose energy density is specified by the scale factor. The scale factor 
a is then the only canonical variable and the Lagrangian for FLRW models can be written in the following way 

m, Hi m [13, ii 


L = — —T7 \a^a — fca] — 27r^a^p, (3-3) 

4Ct 

where p is the total energy density and the dot represents derivative with respect to the cosmic time. Here we are 
neglecting the cosmological constant (which can be incorporated in the dark-energy density). The canonical conjugate 
momentum of the scale factor is 


_ dL Stt 


(3.4) 


and introducing the classical Hamiltonian, H = an a — L, the corresponding constraint 77 = 0 is given by 


- - ^ka + 2n‘^a^p = 0. (3.5) 

Stt a 4G 

This equation can be simplified, yielding precisely the Friedmann equation m- As for the quantum model, depending 
on the factor ordering, we will get different equations. We will verify if the DW condition [13, i-e the vanishing of 
the wave function when the scale factor goes to infinity for a p given by Eq. (1^ . can be satisfied and if it depends 
on the factor ordering choice. 


1. First quantisation procedure (a'H{a,na)'4>{a) = Oj 

The simplest way to obtain the WDW equation with the scale factor as the single degree of freedom is to multiply 
the Hamiltonian constraint, Eq. (13.5L by a and make the substitution tt^ —>■ ffa = The corresponding WDW 

equation is 


, G (P 
3n da? 


V{a) 


tjj{a) = 0, 


(3.6) 


where 


V{a) 


+ 27r^a'‘p(a). 

4G 


(3.7) 


For the case of a DE fluid obeying the equation of state (1231), the wave function of the universe obeys the following 
differential equation 


,9 Cr d 37r 9 ^ A 

37r da^ 4G 


A + A\n( — 
\ao 


ipia) = 0, 


(3.8) 


and for fc = 0, it reduces to 


J d^ 6n^ 4 


A -I- Ain 


Uq 


ip{a) = 0. 


We perform the following changes to achieve a dimensionless WDW equation, 

a ^ 8nG , ^ 8nG , nanHo 


(3.9) 


(3.10) 


"0 


Therefore, the equation (|3.9p becomes 


^ + [Da + ln(u)] ^(u) = 0. 


(3.11) 












Using the WKB (Wentzel, Kramers, Brillouin) approximation method at first order, which is enough to ensure that 
the DW condition is satisfied; i.e. 'I'(m) vanishes for tt —>■ oo, the wave function of the universe can be described by 

IS [Hi 


^(u) [Ua + Ua ln(u)]-’ } , 

where Ci and C 2 are constants and the function Sq{u) can be written according to Eq. (lAllI) as 

i 

S'o(u) Ri -q— [Ua + Ua ln(u)] = . 


(3.12) 


(3.13) 


2. Second quantisation procedure ('H{a,Tta)'4>{o) = 0^ “ Laplace-Beltrami factor ordering 

The Laplace-Beitrami factor ordering is perhaps better motivated for two or higher dimensional WDW equations; 
i.e. with two or more physical degrees of freedom, since it is based on a covariant generalisation of the Laplacian 
operator in minisuperspace [S ] . The starting point for the Laplace-Beltrami factor ordering for the case of a perfect 
fluid specified by the scale factor is to consider the wave equation obtained from the classical Hamiltonian constraint, 
Eq. (13.5L without multiplying by a (c.f. also, for example, Ref. 0). Therefore, 


37r a 


+ 27r^a^/9 


V'(a) = 0, 


here we are already restricting the analysis to a flat geometry (fc = 0). Introducing the change of variable 


(3.14) 


such that 


-2^ = -ff 
a 


'-id' 


'-id' 



“ "da 


9h^ 

4 Oq dx^ ’ 


the expression (13.141) can be rewritten as 


dx'^ 


+ Vix) 


if{a) = 0. 


Using the definitions of the parameters realised in Eq. (I3.10L the effective potential V{x) is the following 

V (x) = q^x^ Ua + ^Ua ln(a;) 

The general solution in the WKB, first order approximation, in this case also verifies the DW condition 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


ipix) « — 
y/V 


Da + -Da ln(x) 




where Ci and C 2 are constants and the function So{x) can be written according to Eq. (lAllI) as 

2 

Soix) R! ?7y [Da + Da ln(a:)]" . 


(3.19) 


(3.20) 


IV. QUANTUM STUDY WITH A PHANTOM SCALAR FIELD 


In this section, we study the quantum behaviour of the LSBR using an approach based on considering a phantom 
scalar field. We consider the Born-Oppenheimer (BO) approximation for the the WDW equation and, as in the 
previous section, we check the fulfilment of the DW criterium as an indication of singularity avoidance. 
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A. Wheeler—DeWitt equation with a phantom field - Review 


Assuming FLRW symmetry the gravitational action with a minimally coupled phantom scalar field (for vanishing 
cosmological constant) can be written as (c.f. Ref. 0): 


S = ^ J [—a^a + ka] dt — J ( 


+ 2V{P) 


dt. 


(4.1) 


The corresponding Wheeler-DeWitt equation 'Hip{a, <f>) = 0, for flat spatial geometry is 


2 6 da^ 


d(jP 


'ipia, p) + a);e®“R {(j))ip{a, </>) = 0, 


(4.2) 


where a = ln(a/ao), = SttG, V{(j>) is the potential for the scalar field, which in our model is quadratic, Eq. (I2.10L 
and we have taken k = 0. Here we have used the Laplace-Beltrami factor-ordering since, as discussed in section III, 
in two dimensions it leads to the Laplace-Beltrami operator which is the covariant generalisation of the Laplacian 
operator in minisuperspace [H^. Apart from this formal motivation, we restrict to the consideration of this factor 
ordering because we have seen on the previous section that a change on the factor ordering seems to leave unaffected 
the conclusion about LSBR avoidance. 

Before proceeding to solve the WDW equation some comments are in order. In the first place, one should keep 
in mind that throughout this paper the validity of the quantum description based on the WDW equation is limited 
to a subset of the minisuperspace corresponding to the validity of the DE dominance assumption (ptot ~ Pde = P 
in Eq. (EU). This implies that we are considering through the paper the existence of a finite minimum for a, 
ttmin > —A/A, as can be obtained by considering the positivity of the energy density given by Eq. (12.41) . In the 
second place, a and p are not independent classically, according to Eq. (12.81) . Naively, this could seem to imply 
that to check the fulfilment of the DW criterium we should search for solutions ipia, p) that decay along that line 
(j) = (a), —>■ 0, when the scale factor goes to infinity. Nevertheless, when the energy density 

approaches a certain quantum gravity threshold, the quantum effects should become important breaking the classical 
constraint p = (a). Indeed, although the potential is classically obtained in section B, a and p are independent 

in the quantum description. Therefore, we must investigate the evolution of the wave function </>) for a —>■ oo and 
arbitrary values of (p. Finally, for a compatible link between the quantum solutions and the classical behaviour, far 
from the quantum effects dominance, the wave function should be described by wave-packets centered in each point 
of the classical trajectory (p = (/i'^*“®®(a) [s^. We will not consider this issue in the present work. 


B. Born-Oppenheimer approximation 

To solve the WDW equation (14.21) we will now adopt the BO approximation which considers that the wave function 
can be factorised into two parts, corresponding to the geometric and matter parts. This is 

'ipk{a,4>) = ipkia,4>)Ckia). (4.3) 

The matter part (p{a, cp) is assumed to satisfy the hypothesis of adiabatic dependence with the scale factor, such that 
the terms containing its derivatives with respect to a can be neglected. The consideration of the BO ansatz (14.31) in 
the WDW equation (14.21) leads to two differential equations, one for the geometric part of the wave function and the 
other for the matter part. Those are 

+ [oy'-VW - E,(a)] w(a,« = 0, (4.4) 


d’^Ckja) 
6 da^ 


+ k'^{a)Ck{a) = 0, 


(4.5) 


where k'^{a) = -^Eki^a). In this approximation, the matter part of the wave function has no backreaction in the 
geometric part. The two equations become decoupled and the matter part has only an indirect influence on the 
geometric part through the eigenvalues Ekia). The validity of this Born-Oppenheimer approximation is discussed in 
appendix C. 
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1. Solving the matter part 


Substituting the potential given by Eq. (I2.10L in the matter part of the WDW equation, Eq. (14.4L we obtain 


d'^ipkia,4>) 

2 dcj)^ 


+ S “nC 


6^6a 


A 


27rGA((/>-</)iy 


- Ekia) > ipkia,(j)) = 0 . 


(4.6) 


This equation is not analogous to the quantum harmonic oscillator but is similar to a repulsor instead. This is a 
consequence of considering a phantom scalar field and rather than a usual (quintessence) scalar field. Therefore, the 
spectrum is continuous in this case and we can consider that Ek{a) = Ek and k‘^{a) = are continuous parameters 
independent of a, where k is not restricted to real values. The solutions of Eq. (14.61) are (see appendix B for details) 

t^i'’"^(a,<^) = TT(/l,±z), (4.7) 


^ (a, 0) = K-^/^W (/3, z) + i(3, -z ), 

and 

(a, </.) = K-^/^W (/3, z) - iK^I^W (/3, -z), 
where W (/3, z) is a parabolic cylinder function [g^, K = vT+'e^^— 


2;i(7rG)i/2 [ 6 




Ek 




and 


z = 


2a^/"e3“/2(7rG7l)i/4 


(</>- <(>i)- 


Thus, the matter part of the wave function is a combinations of these solutions 

ipk{.a,4)) = ciip^^\a,(y)+C2‘ff\a,(iy+C3'^f\a,4y + CA'-p‘t\a,4>), 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 


with q’s constants. As we are dealing with a second order differential equation only at most two of the previous 
functions are linearly independent. 

As shown in appendix B, in the limit a —>■ c» these parabolic cylinder functions go as 

W (/3, z) - cos(e3“), (4.13) 

and 

W (/3, -z) ~ sin(e3“), (4.14) 


Therefore, the matter part of the wave function tend to zero when a —^ c» for arbitrary Ci’s. This imply that the DW 
condition is satisfied if the geometric part of the wave function stays hnite in this limit. 


2. Solutions to the geometric part 
The solutions of equation (14.51) are oscillatory for Ek > 0, 

. y/T2EZ . ^/12 eZ 

Gfe(a) = aie* '‘K “ + 026 “® -'R “ 
and have an exponential behaviour for Ek < 0 

Gfc(a) = 5ie “+626 ‘ 


(4.15) 


(4.16) 


where ai,a 2 ,&i ,&2 are constants. Therefore, to make sure that DW criterion for singularity avoidance is fulfilled we 
impose bi = 0 . 
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V. CONCLUSIONS AND OUTLOOK 

The main purpose of the present work was to analyse the behaviour of the Universe when approaching a LSBR [l[ 
within the framework of quantum cosmology [13, HI] • We have obtained the WDW equation in the context of quantum 
geometrodynamics for two cases: i) a perfect fluid with the equation of state (12. 3|) and ii) a dominant scalar field of 
phantom character whose potential is given in Eq. (12.101) . We applied the DW criterium for singularity avoidance, 
namely the condition that the wave function should vanish in the region of the classical singularity. In the first case, 
the DW condition is satisfied, independently of the two factor ordering considered. This can easily be verified by the 
solutions obtained in the WKB approximation, Eqs. (13.121) and (13.191) . Eor two degrees of freedom, we have also 
shown that it is possible to find solutions obeying the DW condition for any fixed value of the phantom field </> or 
even along the classical trajectory, c.f. Eqs. (HTl) . (ITO . and (HU). These results might point to a possible resolution 
of the LSBR in the quantum description. 

Eor the scalar field case, we used a BO approximation leading to two separate equations, one for the geometric part 
of the wave function and the other for the matter part. In this approximation, the matter part of the wave function 
has no backreaction in the geometric part. The two equations become decoupled and the matter part has only an 
indirect influence on the geometric part through the eigen values E^. The solutions for the geometric part can be 
divergent (which we disregard) and convergent exponentials or oscillatory for Ek < 0 or Ek > 0, respectively. Due 
to the phantom character of the scalar field, the equation for the matter part is analogous to a quantum harmonic 
repulsor. The solutions of this equation are Parabolic Cylinder functions |68|. As it is shown in appendix B in detail, 
these solutions vanish in the limit of arbitrarily large scale factor. Therefore, the DW condition is asymptotically 
satisfied. 

The fundamental starting point for the classical model with the LSBR is the equation of state (lOl) . slightly different 
from that of a cosmological constant. For the case of a perfect fluid, this relation, which is implicit in the energy 
density, is directly inserted in the WDW equation in the matter component. Since in this case the matter component 
has no implicit degrees of freedom, one can say that the traces of the LSBR equation of state are fully encapsulated in 
the geometric part of the WDW. On the other hand, for the phantom scalar field case, the WDW is two-dimensional 
and the LSBR signature can be inserted on the phantom part via the scalar field potential. 

As mentioned in the introduction, the issue of boundary or initial conditions for the wave function is of extreme 
relevance in quantum cosmology and it still is in some sense an open question (s^-lssj]. The no-boundary wave function, 
the Tunnellin g w ave functions, were essentially motivated by the application of quantum cosmology to the very early 
universe (see [56l| and references therein). The no-boundary proposal was in some sense unique since it tried to bring 
together the dynamics and the problem of the initial state. Nevertheless, so far there is no unified theory of boundary 
conditions and dynamics from which that relation can be derived . A more general suggestion, the one we used in 
this work, was first given by DeWitt, in what has been named the DeWitt condition. It is a general condition imposed 
on the wave function saying that it should be zero in the regions corresponding to the classical singularities. Again, 
at the present moment there is no fundamental proof based on quantum gravity that the solutions to the WDW 
equation should obey this condition . If this is found then one could say that the theory excludes the (geometric) 
singularity problems without any ambiguity. 

For genuinely quantum regions, the correspondence between the classical and quantum predictions should in prin¬ 
ciple break down and these are expected to occur at the vicinity of classical singularities. Now, as pointed in the 
introduction, one way for this to happen is to get wave packets solutions (follo wing the classical trajectory) that 
become smeared out in the classical singularity region, effectively avoiding it [441154| and another way is to have a 
decaying wave function in those regions, satisfying the DW condition [s^j- The DW condition, taken as a criterium 
to avoid the classical singularities is based on the existence of square integral functions, and therefore on a consistent 
probability interpretation for the wave function. Accordingly the probability amplitudes for wave packets should 
vanish in the limit when the variables (of the representation) go to infinity. The problem is the fact that these square 
integral functions require an appropriate Hilbert space and it is not obvious that this can always be done in the 
quantum cosmology based on the WDW equation (c.f. for example the references [bOl - f^ for cases where this is 
doable). 

The LSBR is an event happening in the limit when the scale factor goes to infinity, leading to the divergence of 
the Ricci parameter, the Hubble parameter and of course the dominant dark energy density. Can we avoid, with the 
quantum treatment, the region of configuration space corresponding to the divergent scale factor? As mentioned, in 
this work we tried to impose the DW condition on the solutions to the WDW equation in order to achieve this goal. 
What about the other quantities mentioned: how to compute them in the quantum geometrodynamics approach? In 
the quantum description, cosmological quantities such as the scale factor, the Hubble parameter (and its derivatives), 
or the Ricci curvature are represented by operators. Therefore, in principle if we have a normalised wave function 
in minisuperspace we can evaluate the expectation values and probability densities for these observables. Again, we 
emphasise the fact that the interpretation of the results for the expectation values in the quantum cosmology based on 
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the WDW equation depends crucially on a consistent probability interpretation for the wave function and therefore 
on a minisuperspace with a proper Hilbert space nature, which remains an open question. Strictly speaking, the 
DW condition on the wave function is not sufficient to guarantee that all the relevant quantities which diverge in the 
classical model become finite in the quantum description. One should compute the corresponding expectation values 
and/or probability amplitudes to have a complete analysis. In this work, we simply checked if the DW condition could 
be verified. The calculation of expectation values and probability amplitudes requires further investigation since it 
is intimately linked to various open questions in quantum cosmology regarding a Hilbert space structure and the 
classical-quantum correspondence. 

Nevertheless, it can be shown for the case with the scale factor as the single degree of freedom, for solutions Eos. (13.121) 
and (|3.19l) obtained in the WKB method, that the energy density expectation value is finite. If the effective Friedmann 
equation, obtained by replacing the classical quantities by the corresponding expectation values, remains valid, it would 
then imply that the Hubble rate expectation value is also finite in the quantum description. Nevertheless, this is not 
completely clear and requires further investigation, first of all because the calculation of the expectation value of the 
Hubble parameter depends on the factor ordering and second, because this effective Friedmann equation might be only 
valid for the classical regime. In fact, the standard quantum cosmology implied here, based on the WDW equation, 
is strictly speaking a canonical quantisation of the classical theory of GR. In the genuine quantum regions mentioned 
the classical gravitational theory of GR is expected to break and we might need a different theory of quantum gravity, 
although as it was mentioned in the introduction, many results regarding the WDW analysis of late time cosmology 


are expected to hold [ij. Does the (effective) Friedman equation < i/) H^(a, TTa) -0 >= < -0 \p{a)\ 0 >, using the 


solution to the WDW equation, correspond exactly to the classical one? If not, is there any region in minisuperspace 
in which it is valid, in accordance with the correspondence principle? Notice that the above expression contrasts 
with the semi-classical approach in which no quantisation is done in the geometric side of the Einstein equations, 

R/j-v — < 'P T/jy 'P > (c.f. for example Ref. 0). Historically this approach has been comparatively 


more explored, in which the energy momentum is replaced by the expectation value computed from the quantum 
state describing the matter, where W is the wave function representing the degrees of freedom of some quantum field 
theory in the fixed background of curved space-time. These issues regarding the calculation of cosmological quantities 
in the quantum description and their relation to the classical model require further investigation. 

Regarding the topic of future cosmic singularities, there are many open questions. First of all, what kind of 
observational predictions could distinguish between the different models with late-time acceleration and dark energy 
related singularities? Another interesting aspect implicit at least in some of these models is the possibility of matter- 
energy transgressing most or even all of the classical energy conditions. This lead to the idea that these conditions 
can be refined through semi-classical or purely quantum arguments A new door on the physics of phase 

transitions might be open if these non-linear quantum energy conditions (and their interpretations) are profoundly 
explored. 
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Appendix A: The WKB approximation 


1. general description 


For the following second order differential equation 


I?+ ’'<*') 


0(?/) = 0, 


(Al) 
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where 

V{y) = [12a + 61n(?/)], 

the lowest order in the WKB approximation gives [5^ 

V'(y) « 

where Bi and B 2 are constants and 


So{y) = / \/V{y)dy. 


(A2) 

(A3) 

(A4) 


'yi 


The variable y can be written as a positive power of a, according to ()3.10p and (13.151) . Therefore, when a —>■ 00 , 
y ^ 00 . Since V{y) is always positive, the solutions for the zeroth order WKB approximation are purely oscillatory. 
The lower order WKB approximation do not fulfill the DW condition, so we go to next order. 

The WKB approximation in the first order will allow us to have an approximate evaluation of the asymptotic 
behavior of the independent solutions to our wave equation, the general expression is 


^{y) Ki [-V{y)Y 


Bie 


iSo{y) 


Bne 


-iSoiy) 


(A5) 


As we can see, for large values of y, the function V{y) diverges and therefore, the wave-function, tpiy) vanishes. The 
DW condition is satisfied and the Little Sibling event is effectively avoided in the quantum description with one degree 
of freedom, independently of the two factor orderings considered here. Only the parameters s, b and ?) depend on the 
factor ordering: 

1. For the first quantization procedure used in the present work, we use the operator a'H[a,'ka). Then, the 
parameters above are written as follows 

n i. ~ SnalHo 

S = 2, b = nA, y = . (A 6 ) 

2. For the second factor ordering choice, which corresponds to the operator 'H{a,TTa), the parameters are written 
as 


= 1 , 5 : 


2 . TragiLo 


The integral in Eq. (|X4l) is given by 


So{y) = / TO® [Oa -f bln{y)] = dy. 


after performing an integral by parts we get 


So{y) = -k 61n(y)]^ 

s -I- 1 




by^ 


'VI 


' 2 (s + l) 


[OA-k61n(y)] ^ dy. 


(A7) 


(AS) 


(A9) 


The integral of the second term can be done defining an auxiliary variable t = [Oa + &lii(y)]. This leads to 




So{y) = [^2a + dln(y)] ^ - fj 


y/Ybe ^ (®+i) 
2 (s-k 1 )^ 


erfi 


yf+^ 1 

[f^A -k61n(yi)]- 

S ~r -L 


77 - 




erfi 


s -I- 1 
b 

s + 1 


a/Da + &ln(y) 


i/riA -I- 61n(yi) 


(AlO) 


2 (s + l)= t 

where yi is an arbitrary integration constant such that oi (corresponding to yi) is larger than oq. The function 
erh( 2 :) is the “imaginary error function”, which is related with the error function (defined by Eq. 7.1.1 of page 297 in 
Ref. [ 63 ) through erfi(z) = —ierf(zz). The second term on the rhs in Eq. (lAlOl) is much smaller than the first term 
for large values of y (as it can be noted considering Eq. 7.1.23 of page 298 in Ref. [g^), and the third and four terms 
can, of course, be dismissed in this limit. Therefore, we can neglect it for big values of y. 

iiS+1 , 


So{y) ~ [^A + bMy)Y 

s + l 


(All) 


This approximation (lAllI) is then applied in equations (13.131) and (13.191) . where the value of the parameters are 
written in (IMl) and (TATI) . for the first and the second quantisation procedures, respectively. 
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2. Validity of the WKB approximation. 

The wave equations we have obtained in Eqs. (13.191) and (I3.12|) using the two different factor-orderings are both of 
the form 


[dl + fi^giy)] ilj{y) = 0, (A12) 

for equations of the form (IA12I) . the validity of the WKB method in the zeroth order approximation; i.e. the solutions 
(13.121) and (|3.19l) without the global multiplicative prefactors, is ensured if 


1 

V 


1 dg{y) 

2g{y)^ dy 


< 1 . 


(A13) 


On the other hand, for the first order approximation the inequality that must be satisfied is 


1 - 4g(2/)g(y) 

72 16g3(y) 


(AM) 


where dot stands for a derivative respect to y. This condition must be fulfilled for the validity of the equations (13.121) 
and (13.191) . These inequalities are obtained as conditions in order to neglect the remaining terms which are absent in 
the corresponding approximations. Being g{y) = [7 -I- /31n(y)] we get 

9{y) = {/3 + 2 s[ 7 -f ^ln(y)]}y2®"\ g{y) = {( 2 s - 1 ) {/3 -f 2s [ 7 -f/?ln(y)]} 2s/3} (A15) 

so for the zeroth order approximation we have 


1 

77 


{/3-f2s[7-f/31n(y)]} ^_(,^i) 
2 [ 7 -f/31n(y)] = 


< 1 , 


and for the first order approximation we get 


y-2(s-tl) 

b + PHy)? 


5 {l3 + 2s[y + ^\n{y)]f 
16 [7-f/31n(y)] 


i{(2s- l){/3-f 2s[7-f/31n(y)]}-h2s/3} 


< 1 


(A16) 


(A17) 


which is clearly verified when y —> 00 because —1 < s, for both cases and for the two quantization methods. Here 7 
and /3 are constants which for the zeroth and first orders are found in Eqs. (13.111) and (I3.18p . respectively. 


Appendix B: Detailed calculations for the matter part of the WDW equation 

The matter part of the WDW equation, given by Eq. gH), can be rewritten as 

d^ipk{a,z) , fl 2 


dz^ 


where 


and 


/3 = - 


+ [-Z - /3 ipk{a,z) = 0, 


Ek 


2/i(7rG)^/2 6 

2ao/^e3“/2(7rGA)i/4 








(Bl) 


(B2) 


(B3) 


Eq. (IBII) corresponds to Eq. 19.1.3 of page 686 in Ref. [ 6 ^. Thus, its solutions are the following parabolic cylinder 
functions: 

= W {13, z), = W {P,-z), = E{P, z), and = E* {P, z), (B4) 
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with 


and 


E{a, (f>) = [P, z) + iK^/'^W (/3, -z ), 


(B5) 


K = v'l + 


(B6) 


We are interested in the behaviour of these functions for arbitrary values of the field and very large values of a. It 
can be noted that /3 —>■ — oo when a —> oo. Defining 


X = - 4/3 




l2TTGA{(j)-PiY 


A 


Ek 


T 6 ^ 6 a ’ 


i-QC 


(B7) 


it can be noted that X —>■ oo when a —>■ oo. Thus, we focus our attention on the case /? < 0 and 0. The real 

solutions for this regime are given by (Eqs. 19.23.10 and 19.23.11 of page 694 in reference [ 13 ) 


and 


with 


W (P, z) = \/2Ke^'' cos(7r/4 + 0 + Vi), 
W {P, -z) = \/2lKe'’'' sin(7r/4 + 0 + Vi), 


e = ]-zx-p\i,(^^^ 

4 V^/W\) 

^3^3a 


, s/7ga{4, - mA.ga(4. + 


A Ek 

In 

V TT GA{(j) — 4>i) + ^ 2 ttGA{ 4) — piY + 

6 OgC®" 

/ ^ 



Y 6 aQe®“ 


(B8) 

(B9) 


(BIO) 


and Vi and Vr are defined by series that go as Vr « —1/2 InX and Ui ~ 0 for very large values of X (see Eq. 19.23.4 
of Ref. [63). Moreover, using Eqs. (IB6I) and (IBIOI) . it can be seen that 6 diverges as e^“ and R —>■ 1 for a —>■ oo. For 
large values of a, we can approximate Eq. (IB8I) as 


Thus, we have 


and 


W{P, z) ft! y/2X-^/'^cos 


TT 1 

- + -zX-p\n. ^ 
4 4 


/ z + x ' 


W(P, z) ~ e-3“/^cos(e3“) ^ 0, 


W(P, -z) ~ e-3“/Sin(e3“) ^ 0, 


(Bll) 


(B12) 


(B13) 


for a —>• oo. On the other hand, taking into account Eqs. (IB12I) and (IB13I) into Eq. (IB5I) . one can conclude that the 
imaginary solutions also vanish in the asymptotic limit. 


Appendix C: Validity of the Born—Oppenheimer approximation 

Consider the WDW equation (14.21) for the phantom case, 


2 6 



da^ 


+ ‘fk 


da^ 


dGkd(pk \ 
da da J 


+ + ay‘^V(P)Gk‘Pk = 0 . 


(Cl) 
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The BO approximation assumes that the first and third terms of this equation can be neglected as compared with 
the second term. That is, it considers that the variations with respect to geometry of the geometric part, Ck^o:), are 
much more important than those of the matter part, and that, therefore, there is a negligible backreaction 

of the matter part on the gravitational part. 

The validity of this approximation can be explored using expansions with respect to k (see Ref. and references 
cited there). For the geometric part, from Eq. (I4.5|l . we get 

(C2) 

with ' = djda. The matter part is given in terms of parabolic cylinder functions. In order to get its expansion in 
terms of k and A in a simple way, we use the approximation for large values of the scale factor given by Eq. (IBllj) . 
We should also keep track of the dependence on a. Noting that 

/3 - 0 - 0{K-^A^/^al) e^°‘, z - X ~ , (C3) 


we get 


ipkia, 4>) ~ E{a, (j)) = V 2 exp [h{a, 4>)] 


where 


h{a) ~ O (1) a + O (^K e^“ 


and the quantity k 


is dimensionless. Thus, we have 


Ck'^k 

Ck<fk 


OiK^) 


h^ + h 


OiK^)+0{KA^/^ a3)e3“ +C>(Aa[;)e®“, 


(C4) 


(C5) 


(C6) 


and 


Ck^k 

Ck^k 


0(k) h - 0{k^) + OiA^^^al) 6^“. 


(C7) 


Thus, the B.O. approximation is fulfilled as long as AaQe®“ <C i.e. this approach is valid for sufficiently small 
value of A, i.e. the fluid (12.31) is close enough to a cosmological constant and for large enough values of a but not 
infinite value of a, that is within a semiclassical regime. 
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